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ABSTRACT 


Radar imaging is an area of tremendous interest as radar-based systems are 
perhaps the only all-weather, long range remote sensing systems. However, radar’s 
continued utility and application in wide-ranging areas is fundamentally dependent on the 
ability to produce high-quality, artifact-free imagery. To date, the use of radar to identify 
and image moving objects remains of great interest, and it is well known that motion in 
the scene gives rise to mispositioning or streaking when target motion is not addressed. 
Many techniques have been developed to handle moving objects, but these techniques 
make use of the start-stop approximation, in which a target in motion is assumed to be 
momentarily stationary while it is being interrogated by a radar pulse. 

A new linearized imaging theory that combines the spatial, temporal and spectral 
aspects of scattered waves has been developed. This thesis studies the performance of 
this new imaging scheme via analysis of the point spread function. It is shown that the 
imaging PSF localizes, and is translation invariant in phase-space. It is also shown that 
the behavior of the imaging system is dependent on the aperture geometry and choice of 
radar waveforms transmitted. 
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I. 


INTRODUCTION 


A. RADAR IMAGING 

Sensor systems that can detect, locate and identify targets at great distances and in 
all kinds of weather have well-recognized utility. Radar-based systems are capable of 
performing such functions over long-range, all-weather and all-day (day or night) unlike 
many optical counterparts. This is because the wavelength of radar signals makes them 
relatively unaffected by atmospheric and weather-induced attenuation. Herein also lies 
the problem: image resolution depends upon signal wavelength. The resultant side effect 
is that resolution suffers and radar-based imaging will be fundamentally less effective 
compared to other systems based upon relatively shorter signal wavelengths. 
Nevertheless, radar waves scatter strongly from objects whose size is on the same order 
as the wavelength. This means that radar is sensitive to objects whose length scales range 
from centimeters to meters. Many objects of interest are in this range. 

Many research and development efforts have been expended to improve the use of 
radar for imaging. It has been reported that resolution on the order of centimeters can be 
achieved with spotlight-mode operation [15]. While the developed techniques have been 
instrumental in making radar imaging a more viable application, the parallel advent of 
fast digital computers and leaps in computational prowess, have enabled practical 
systems to be fielded. This has generated great activity in areas such as interferometric 
synthetic aperture radar (SAR) and bistatic SAR, due to interest in operational concepts. 
The key to success in all of these applications is the ability to produce high-quality, 
artifact-free imagery. 

B. IMAGE ARTIFACTS 

The subject of radar (and radar imaging) has evolved considerably since the 
1950s. Consequently, many implicit assumptions are embedded within current imaging 
techniques. These simplifications or system limitations are of interest as they are 
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frequently also responsible for the image artifacts discussed in the literature. The start- 
stop approximation is central to the thesis and holds the fundamental difference between 
the new imaging scheme developed in [2] and current techniques. 

1. Scattering Model - Born Approximation 

Historically, scattering models have been built upon a foundation which views 
targets as being composed of simple, non-interacting point scatterers. When the imaging 
scene consists of multiple targets, accounting for multiple scattering (in which waves 
bounce back and forth between the targets) can become extremely complex. For this 
reason, the ‘weak scatterer’ or Bom approximation is typically invoked. A useful 
consequence is that the imaging problem becomes linear and more tractable. 

However, the Born approximation may result in image artifacts as actual 
“multiple-scattered” waves are received later than expected from the single-scattering 
model. These waves are then interpreted as having come from a more distant target 
structure. In addition, increased capabilities in radar resolution and sensitivity and 
parallel advancements in computational technology have exposed the inadequacies of 
such simple scattering models. With the ability of radar systems to resolve and image the 
components (subscatterers) that lie within the support of a traditional target, it is 
inevitable that new radar target models that can describe these multiple scattering effects 
accurately (whether multiple targets or multiple components of a single target) are 
required. Current approaches typically utilize correction techniques to enhance imaging 
algorithms based on the low order ‘weak scatterer’ model. 

2. Measurement Systems 

There are inherent limitations in any data that are measured by radar systems. 
Such data are usually corrupted by noise and are always limited by bandwidth. 
Consequently, inducing unwanted image artifacts and reduced image resolution. 
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3. 


“Start-Stop” Approximation 


The radar uncertainty principle bounds the accuracy with which time delay and 
frequency shift can be simultaneously estimated using a single pulse. This limitation may 
be mitigated by the use of follow-on pulses to further probe target characteristics. Most 
modem radars use a train of high range-resolution waveforms to enable pulse-to-pulse 
velocity estimation. Such waveforms allow the use of the start-stop approximation, 
which assumes that the target is stationary during the time the radar pulse illuminates the 
target, and moves only between pulses. When individual pulses are short, and targets of 
interest are small and move slowly compared to the speed of light, this is almost always a 
good approximation. Under these circumstances, good velocity estimates are obtained 
from the range-rate (two accurate range measurements, separated by a time interval). 
Unfortunately, no Doppler information is actually measured by the systems. Range and 
range-rate information are then used to form spatial images. 

In order to achieve the necessary signal-to-noise ratio (by coherent integration), 
long pulse trains are required. Consequently, during the coherent processing interval, the 
stationary target assumption may no longer hold. Scatterers may move to different 
positions, and when they move beyond a resolution cell, the result is a blurred image. 

C. MOTIVATION - IMAGING MOVING TARGETS 

In recent years, a number of attempts to develop imaging techniques that can 
handle moving objects have been proposed. Space-time adaptive processing (STAP) is a 
signal processing technique that was originally developed for detecting slow-moving 
targets using airborne radars. STAP uses multiple-element antenna array, together with 
real-aperture imaging techniques to produce ground moving target indicator (GMTI) 
images [23]. This method described in [24] combines SAR (designed to image stationary 
scenes) and GMTI processing to develop new techniques for detecting slow-moving 
surface targets that exhibit start-stop maneuvers. The method described in [11] discusses 
the extension to non-sideways looking array radar. 
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The method described in [10] describes the concept of velocity synthetic aperture 
radar (VSAR). VSAR is a multi-element SAR system involving conventional processing 
to form an image at each element; except the image phases are preserved and compared 
to deduce target velocity. However, VSAR processing assumes that the scatterers remain 
in a single resolution cell throughout the integration time. 

Time-frequency methods for signal analysis allow decomposition of the 
frequency spectrum of time-varying signals in shorter time-frames, and provide an added 
dimension for examining the dynamic behavior of the signal as it varies over time. The 
size and shape of the weighting window function can be altered to fit the specific analysis 
requirement for the signal. In recent years, further extension of this time-frequency 
method to ISAR has been advocated in the area of image generation, motion 
compensation and micro-Doppler target vibration studies [Refs. 8, 20 and 21]. Chen [20] 
has reportedly used the time-frequency transform (TFT) to construct inverse synthetic 
aperture radar (ISAR) images of targets with high rotation rate by helping to overcome 
the migration of individual scattering points from one range cell to another. 

However, all of these techniques rely on making the approximation that a target in 
motion is assumed to be momentarily stationary while interrogated by a radar pulse. It is 
apparent that the utility of an imaging approach that can accommodate target motion 
during the imaging process. This constitutes the focus of the thesis. 

D. OBJECTIVE 

The objective of this thesis is to study a linearized imaging theory developed by 
Cheney , Borden [2] for an imaging region containing moving objects. The physics 
involved, and the approach for addressing image artifacts common to targets moving in 
an unknown fashion will be discussed. Analysis of the corresponding point-spread 
function will be carried out to examine the localization in phase-space (position-velocity). 
The behavior (position, and velocity resolution) of the imaging system will also be 
examined under different aperture geometries and transmitted waveforms. 
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E. THESIS ORGANIZATION 


The thesis is organized in the following manner: 

Chapter II serves to develop the appropriate radar imaging theory by looking at 
the basic building blocks necessary for imaging - radar operating principle, scattering 
models, correlation reception - and how these components translate to imagery. Finally, 
the thesis includes an overview of the inverse nature of the radar imaging problem to 
provide better appreciation of issues pertaining to practical measurement systems and 
how information is extracted from measurements to obtain the best estimates of the target 
object. 

Chapter III applies the concepts put forth in Chapter II to the imaging scenario. 
The discussion here is based on the work of Cheney, Borden [2] in imaging moving 
targets from scattered waves. The effect of multiple moving targets on the scattering 
model is described. The approximations invoked, specifically slow-moving, narrow-band 
and far-field, will be examined. A scattering model to be used in the simulation program 
will be presented. 

Chapter IV discusses the radar ambiguity function, and the relationship between 
estimation errors and the transmitted waveform. Particular emphasis will be given to 
image quality analysis and the derivation of the point-spread function to quantify the 
performance of the imaging scheme developed in Chapter III. 

Chapter V examines further implications of the transmitted waveform through the 
pulse compression perspective. In particular, the Chirp signal is discussed. 

To demonstrate the behavior of the imaging scheme, a simple simulation program 
is developed for modeling the point-spread function. Chapter VI discusses the results of 
the simulation, specifically the effect of different aperture geometries and transmitted 
waveform characteristics. It also concludes and summarizes the major findings of the 
thesis work and provides recommendations for future work. 
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II. IMAGING THEORY 


Prior to deriving the radar data model for imaging moving targets, Chapter II 
examines the formation of spatial images. It begins by examining the physics behind 
extracting range and range-rate (Doppler shift) information by a radar system. A simple 
range-Doppler imaging example is presented. The chapter concludes with a brief look at 
the inverse nature of the imaging problem and characteristics of measurement systems, 
and from there some of the issues/challenges pertaining to radar imaging. 

A. BASIC RADAR PRINCIPLES 

Most fielded radar systems are pulsed systems, as opposed to stepped-frequency 
systems. Radar systems measure the strength of the backscattered field and the round trip 
delay time of transmitted signal pulse reflected from distant objects. Since the radar 
pulse travels at the speed of light, it is practical to use the measured time delay for the 
round trip to determine the range to the reflecting object. This is the basic principle of 
radar theory and we can determine the range R from R = ct / 2 (where r is the round trip 
time delay). 

Several concepts will be useful in appreciating how radar signals convey 
information for the imaging output. The first is concerned with the kind of signals that a 
radar uses to generate the waveforms it transmits. Radar signals are real-valued functions 
of time but it is convenient to express the waveforms as complex-valued functions. 
Radar systems are able to measure these signals (in both amplitude and phase) as 
functions of time. The imaginary part of the complex signal provides the phase of the 
waveform revealing considerable information about the target. 

Radar waves are electromagnetic waves governed by Maxwell’s equations. The 
waveforms generated from radar signals satisfy the wave equation and the radiation 
condition. These properties ‘code’ the echo (transmitted wave scattered from a radar 
target) waveform with important target specific information such as range, speed and 
bearing. 
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The resolution in the down-range direction is governed by the pulse width, and 
correspondingly its pulse bandwidth. Typical carrier frequency operation bands occur 
within the low attenuation atmospheric windows ~ 1-40 GHz. The frequency band of the 
system, in practice, signifies that radar signals are always of finite duration. Intuitively, a 
higher bandwidth will imply finer resolution for imagery in the range direction. More 
about the concept of “resolution” will be addressed in Chapter IV. 

B. SCATTERING OF ELECTROMAGNETIC WAVES 

Understanding the behavior of electromagnetic waves, specifically propagation 
and scattering of these waves, is essential to appreciating the information that can be 
extracted from the transmitted signals. The next two sections will elaborate this further. 

The discussion of [18] describes a simple one-dimensional scattering model for a 
moving PEC 1 plate in free space (a good approximation for typical radar wave 
propagation). An elementary consequence of Maxwell’s equations is that an 
electromagnetic wave (in the far-field) must have the electric field E\ magnetic field 

A 

and wave vector (direction of wave propagation) k mutually perpendicular. For the 
transmitted wave with field components orthogonal to the x direction (i.e. in the y-z, 
plane), the boundary condition on the PEC results in the following relationship between 
scattered and transmitted electric field: 

fr„g(t + R(t)lc) = -ff°J(,-R(,)lc) 

where the scattered field consists of a left-travelling wave (in the negative x direction). 
The position of the moving plate is given by R(t ). Assuming for simplicity that the plate 
is undergoing linear motion 2 (with velocity v in the x-direction), so that 


' A perfect electrical conductor (PEC) allows the charges to move freely and instantaneously in 
response to a field; consequently the fields inside a PEC are zero. The resulting boundary conditions are; 
the tangential components of the electric field must be zero; and the tangential components of the magnetic 
field are related to currents flowing on the surface of the PEC. 

2 Even if the motion is not linear, the linear expression R(t) = x+vt can be considered to be the first two 
terms in a Taylor expansion that is valid for small t (i.e., for short pulses). 
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R(t) = x+vt 


x=R (at boundary) 


and making use of the substitution u = t + R(t ) / c , we can 


solve for t in terms of u (at the boundary): 

u - R / c 

t = - 

1 + v/c 

and obtain the expression for the scattered field 

&° catt g(t + R(t) / c) = -&IJ [a(t + x/c-R/c)-R/c] (2.1) 


where a is the Doppler scale factor 


1 -vie 

a =- 

1 + v/c 


( 2 . 2 ) 


For a transmitted field consisting of a signal waveform s(t ) mixed with a carrier 
wave cos (at t) , the scattered signal received at the antenna (x = 0) is 

P rec (t) ~ s(t - 2 R / c)c os (oj () [(1 - 2v / c)(t -R/c)-R/cJ) (2.3) 

We have used the fact that s(t ) is slowly varying to approximate the argument 
a(t - R / c) - R / c by t -2R/ c , and also the fact that for a scatterer moving slowly, so 
that vie is small. In addition, we can expand the denominator of Equation (2.2) in a 
geometric series and approximate a by 1-2 vie. The multiplier of 1 in the argument of 
the cosine function, is the carrier frequency (relative to the transmitted carrier frequency) 
shifted by an amount 

2v 

u =-(2-4) 

c 


This is the Doppler shift. As in the stationary scatterer case, the received signal is time- 
delayed, but here it is also a Doppler-shifted version of the transmitted signal. 


C. CORRELATION RECEPTION - RADAR DATA MODEL 


Detection of signals in noise is fundamental to radar operations. A radar detects 
the presence of an echo signal reflected from a target and extracts information about the 
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target for uses such as imaging. The conservation of energy requires that the wave 
intensity from a finite-sized antenna decreases in strength as the inverse square of the 
distance. Consequently, the round trip power of the signal received by the same antenna 
for small targets decays by a factor of 1 / R 4 relative to the signal that was transmitted. 
For typical radar operational distances, the received power can be as low as 10 18 Watts 
consequently. Thermal noise in the receiver presents a problem for target detection, and 
a greater problem for image formation. 

While increasing the power of the transmitter (or pulse integration) can help 
increase the signal-to-noise ratio (SNR), there are practical limits to relying solely on 
these approaches. Key amongst these limits is the decrease in maximum unambiguous 
range, and also the fact that targets move during data collection for pulse integration, 
altering the phase of the scattered field. It should be pointed out that the ensuing 
discussion is premised on the correlation comparison being performed coherently such 
that the phase of the transmission signal is preserved in the reference signal. 

In practice, a solution involves transmitting long coded pulses, together with 
appropriate signal processing techniques called matched filtering or correlation 
reception. This is known as pulse compression (see Chapter V). The correlation receiver 
(or matched filter) is an important example of a radar signal processor for detection of 
desired signals. The output of the matched filter is the cross-correlation function of the 
received signal and the expected received signal (which is related to the transmitted 
signal). Hence, the mathematical equivalence between the correlation receiver and the 
matched filter receiver, and it is possible to implement the matched filter as a correlation 
process. 

A radar obtains information about a target by comparing the received echo signal 
with the signal that was transmitted. Correlation reception presents an intuitive approach 
to develop the standard radar data model used in imaging. The signal processing problem 
is that of optimal detection in additive noise. Assuming Gaussian noise statistics, 

s (t') 

maximum likelihood processing (conditioned on the measured signal rec J ) selects the 
best guess estimate for the radar signal associated with the field reflected from a target 
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s (t^) 

scatt y 'that maximizes the conditional probability. Reference [19] shows that such 


processing (averaged over all time) reduces to finding the function 


,(0 


that minimizes 


f Is (t') 2 dt'+ f 5 (t ') 2 dt'- 

J | rec x 7 J scatt v 7 


2Re f 5 (t')s* (t')dt' 

J rec v ' scatt v x 


(2.5) 


S (t ^ iS 1 (t ^ 

The objective is to determine scan ' y from the random measurements of 7 . 
The first two terms are system-scenario specific. Searching for the best minimizing 
function hinges on the last term and can be made more efficient by restricting the search 
to a few parameters derived from some model. In this case, the natural model is based on 

s it') 

the scattering interaction between the interrogating field and the target. If denotes 

an incident pulse transmitted by the radar, then the linear radar (weak-scatterer) scattering 
model follows by superposition 


s sca«W = [f T Kc( f - T)e iv(t - T) drdu (2.6) 


where p(u, r) is the target reflectivity density for a radar system that measures two 
parameters (range and radial velocity). This reflectivity density is defined in such a way 
that p(u,r)drdu is proportional to the field reflected from the target at range between 
cr/2 and c(r + dr) / 2 with Doppler shift between u and u + do. 

From Equation (2.6), the scattered field is naturally modeled by the two 
parameters z and v . The simplified search space will be parameterized by the time shift 
r and the frequency shift u. Considering the last term of Equation (2.5), and 
maximizing for all z, o in the scattered model, the correlation integral is 

oo 

tj{d, z) = | s rec (t')s* nc (t'- z)e~ M, '~ T) dt' (2.7) 


The correlation integral seeks out the component of received signal that “matches” the 
time-delayed, frequency-shifted version of the transmitted signal. Unlike signals will be 
suppressed consequently; correlation receivers attempt to find the z and v that 
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maximize the real part of Equation (2.7). The expected output of the correlation receivers 
is obtained by substituting Equation (2.6) into the correlation integral of Equation (2.7) 


rj(o,r) = JJJ p(v',T')s jnc (t'- r)s* nc (t'- r')e w(t ' ^e^'^dt'dr'do' 


( 2 . 8 ) 


Recall that the signal-processing problem was detection of signals in noise, hence 
the received data s m (t) consists of both the scattered signal and a random noise 

component. The correlation noise term, which is the correlation integral between random 
noise and the time-delayed, frequency-shifted version of the transmitted signal, has been 
dropped from Equation (2.8). The correlation noise term which measures how well the 
scattered signal matched the random noise will usually be small for coherent systems. 
Rearranging Equation (2.8), 


oo 

Ti(v,T) = || p(u\r')z(o-o',r- 

-oo 


i—( u+v')(t-t') 

r')e 2 dr'do' 


(2.9) 


where 


CO . . 

Z(o, r) = J s jnc (t'~ - T)s* nc (t'+ - z)e wt dt' (2.10) 

-oo 

Equation (2.9) is the standard radar data model and expresses the output of the correlation 
receiver as the convolution of p and z ( U P t0 a phase factor). The function j(t>, r) 
defined by Equation (2.10) is the radar ambiguity function. Radar imaging is closely 
related to the radar ambiguity function that characterizes the accuracy that the target 
position and radial velocity are estimated from the radar data (see Chapter IV). 

D. ONE-DIMENSIONAL (HIGH RANGE RESOLUTION) IMAGING 

One-dimensional imaging involves the generation of range profiles, which can be 

used to articulate target substructure. For high range resolution (HRR) techniques, the 

transmitted pulse’s instantaneous range resolution is smaller than that of the target, and as 

it sweeps across the target it sequentially excites the target’s scattering subelements 

which re-radiate energy back to the receiver. When these scattering subelements are non- 
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interacting and point-like, the scattered pulse will be a sum of damped and blurred 
images of the incident pulse, which are shifted by time delays that are proportional to the 
subelement’s range [1], Figure 1 shows an example of one-dimensional images created 
by HRR radar systems. It can be observed that the profile displays additional peaks 
outside of the target support, clear evidence of the problem associated with a single, non¬ 
interacting scatterer approximation. 

y (meter's) 

o 10 20 30 40 



Figure 1. Example of a range profile from a B-727 jetliner. The top view (with 
orientation at time of measurement) is displayed beneath [1], 
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The signal received from a distribution of non-interacting, stationary targets can 
be modeled as 

oo 

s rec (t)= J p(T')s jnc (t- z')dz'+n(t) (2.11) 

-oo 

where n(t) represents an additive and random noise term. Application of a matched 
filter to obtain the radar data model used for imaging 

oo oo oo 

r/(z) = J s rec (t')s* nc (t'- z)dt' = J s* jnc (t'- z) J p(z')s jnc (t'~ z')dz'dt' + noise term 

—oo —oo —oo 

oo 

= JJ s* nc (t'~ z)s jnc (t'~ z')dt'p(z')dz' + noise term 

—00 

Making the substitution t" = t'-z', and neglecting the noise term (negligible with 
matched filtering) gives 

oo 

7(0 = j Z(r~r')p(z')dz' (2.12) 

-oo 

where % is the autocorrelation 

00 00 

Z(r)= J s*(P- z)s(t ")dt " = J s\t')s{t'+ z)dt' (2.13) 

-00 -00 

Equation (2.12) represents a one-dimensional image and shows how the true distribution 
p is related to the radar data model. The image is the convolution of p with %, which in 
high-range resolution (HRR) imaging is also commonly referred to as the point-spread 
function (PSF). The point spread function obtains its name from the fact that if p consists 

of a single point d(t), then the image is rj(t) = ^ %(t-t')S(t')dt'= %(t). Thus, % 

quantifies the degree to which a single point appears spread out in the HRR image. 

The down-range profile can be affected by a variety of factors such as target 
aspect angle, position of the scatterers or masking of scatterers by other parts of the 
target. Additionally, while the use of short pulses enhances the resolution, it concurrently 
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leads to large bandwidth requirements. Wide bandwidth can increase system complexity 
and increase the likelihood of interference from other emitters in the electromagnetic 
spectrum. A short-pulse waveform also provides less accurate radial velocity 
measurement, a natural consequence from the properties of the radar ambiguity function 
(more details on radar resolution and the ambiguity function in Chapter III). An 
important implementation limitation to highlight for practical radar applications is the 
required high peak power to transmit short pulses over long ranges. High peak power 
transmission can be problematic, especially at high frequencies (due to the small 
waveguide dimensions and the small area of the anodes of the microwave sources). 

Radar target recognition using only range profiles has limited applications. This 
is because a range profile will not be able to distinguish cross-range target structures. All 
scatterers located at the same distance from the radar will reflect energy back to the radar 
with the same time delay. Hence, when the radar illuminates many distinct targets at any 
instant, meaningful interpretation of on-scene target (multiple) disposition cannot be 
derived based on a set (single) of range-only data. 

E. TWO-DIMENSIONAL IMAGING 

For more effective target interpretation, additional information is needed on top of 
the range profile. This information can be in the form of high-resolution cross-range 
profile, Doppler profile, or simply the “triangulation” of range profiles. 

As a simple illustration for extension of the radar imaging concept to two 
dimensions, “triangulation” of different sets of range profiles will be examined (Figures 2 
and 3). This approach enables the determination of cross-range target structure while 
using only HRR radar systems and relies on collecting multiple sets of range profiles 
from different target orientations, processing them, and subsequently synthesizing a two- 
dimensional image. 

Consider a set of three point targets with the radar located at the same distance 
from targets 2 and 3 (see Figure 2). When the geometry is oriented as shown, the return 
echo will only indicate two targets. Hence, ambiguity exists when targets lie along the 
bands of constant range from the radar. 
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Target 1 

o 


Bands of constant 
range from radar 


Target 3 


o 


o 


Target 2 


Single pulse ambiguity 


Figure 2. Ambiguous scenario from a single radar pulse [19]. 


With multiple sets of data from different directions, triangulation will allow for 
gradual buildup of the relative positions of the three targets (Figure 3). The range 
profiles are swept in the cross-range direction to form bands of constant range from the 
radar, which represent the possible locations of the target scatterers. These bands are 
then superimposed and the crossing points are used to determine the scattering center 
locations. The success of such a scheme hinges fundamentally on the ability to correlate 
the various constant range bands. Otherwise, the intersection of the swept lines will be 
erroneous, resulting in image artifacts. 


Multiple pulse “imaging” 


Target 1 

o 


Target 3 


o 

Target 2 


o 



Figure 3. Cross-range information obtained from range profiles [19]. 
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The signal due to a moving target is a time-delayed, Doppler-shifted version of 
the transmitted signal. The measured echo signal from a single point-like moving target 
can be expressed as 

s rec (t) = ps inc (t - z)e~ w(t ~ T] + n(t) (2.14) 

With two unknown parameters to estimate, namely the delay z and the Doppler 
shift o , a single matched filter will no longer suffice. A set of matched filters, one for 
every possible Doppler shift, will have to be applied 

oo oo 

r/(z,u) = J h v (t - t')s rec (t')dt' = pj h o {t- t')s. nc {t'~ z)e~ wU ~ T) dt ' + noise term 


To estimate z and v , find their values where I //1 takes on its maximum. For a 
distribution of non-interacting moving targets, the output from the filter bank is 

oo 

T]{T,u) = J s* nc (t'~ T)e~ w(, '~ T) s rec {t')dt' 

-oo 

oo oo 

= J s* nc (t'— z)e^ w{, ^ x) JJ p(r',u')s jnc (tr')e~ w(r T) dz'do'dt' 

-00 -00 

00 

= JJJ s* nc (t'~ z)e-'° (t '- T) s mc (t'- z')e~ ivXt ~ r) dt'p(z',u')dz'du' (2.15) 

-00 

00 

= JJJ s*Jt"+ z'- z)e- ioa " +T '- T) s inc (t")e- ivxn dt"p(z\v')dz'dv' 

—00 

00 

= JJ X( r ~ t\v- u')e~ w(T ~ T) p(z\u')dz'do' 

-oo 


where 


Z(t,v) = J s* mc (t"+ z)e w(, ) s(t")dt" = ^s* nc (t'+ z)s. nc (t')e lv, ‘dt' (2.16) 

-oo -oo 

is the radar ambiguity function. Similar to the one-dimensional case, the ambiguity 
function can be interpreted as an imaging point-spread function for range-Doppler 
imaging. Analysis of the ambiguity function determines the fidelity of the range-Doppler 
imaging process, and consequently the accuracy of target range and velocity estimation. 
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1 . 


Range-Doppler Spatial Image 


With estimates for range and range-rate from radar returns, this information can 
be interpreted to form spatial images. The following treatment illustrates this, provided 
the relative motion between the antenna and the target is known. Consider a simplified 
two-dimensional (2D) geometry between the radar and target as shown in (Figure 4). 
The target is assumed to have translational and rotational motion only in the two- 
dimensional plane relative to a stationary radar platform. 

Consider a rotating target [18]. For simplification, assume that any overall 
translational motion is removed. 



Figure 4. The geometry for locating a point target known to be rotating [18]. 

Consider a point target rotating counter-clockwise in a geometry in which the 
antenna is in the plane of rotation (Figure 4). The coordinate origin coincides with the 
center of rotation; and the center of rotation is a large distance R from the radar antenna. 
The large range approximation means that the incident wave is approximately a plane 
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wave, consequently, all horizontal points are (approximately) at the same range. If the 
target is positioned at (x,y) = ^(cos^sin#), then its range is 

R = R Q + y = R 0 + a sin 6 (2.17) 

If its angular velocity is §*■= clO I dt, then its velocity vector is a$-sin#,cos#), its 
tangential speed is I v 1= a§*. and the down range (y) component of the velocity is 
$=l v I cos 6* = a$cos6. However, since x = ac os 6*, this down-range velocity can be 

written $= x$*. where x is the cross-range component of the position. Note that for 
x > 0, the down-range velocity is negative, whereas for x < 0, it is positive. 

This means that since a radar system can measure range ( R (] + y ) and range-rate 

($=x$, it can determine the (x, y) location of the point in space (provided 6 is 
known). In particular, the radar measures the time delay t = 2R / c and the Doppler shift 

2$ 2x$ 

v d = - u o = - u o 

c c 

which means that the coordinates of the point are given by 

{x - y) = { - v tb T-v (2A8> 

This shows how the range-rate radar measurements can be used to estimate 
Doppler and cross-range information. Further, the simplification introduced highlights to 
the difficulties in practical implementation of the radar scattering model. Estimates of the 
scattered field requires any variations in range ft to be estimable. While it is reasonable 
to assume (in the one-dimensional imaging methods) ft as constant over the data set, the 

two-dimensional imaging schemes invariably will not. R will generally vary due to the 
target (translational) motion that occurs during measurements made while the target 
rotates. The effect of failing to correctly account for this translational motion will be the 
introduction of phase errors into these data, and a source of image corruption. 
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F. RADAR IMAGING - AN INVERSE PROBLEM 

A discussion of imaging theory will not be complete without addressing the 
inverse nature of the imaging problem and its associated issues. As the thesis is focused 
primarily on developing and analyzing the point spread function rather than imaging, the 
intent of this section is to highlight the key characteristics of such inverse problems - 
issues pertaining to practical measurement systems and how information is typically 
extracted from measurements to obtain the best estimates of the target object. 

The radar can be well-approximated by a linear, shift-invariant system. This is an 
important characteristic, as it states that the object function f(x) at each point x has a 
corresponding output image m(x ). Further, the system will measure objects the same no 
matter when they start (or, equivalently, where they are located). Radar imaging will not 
be meaningful if the delayed version of the target looks different. A general formulation 
of the imaging problem can be achieved by considering an imaging system represented 
by the functional operator (“kernel”) k that describes how the measurement system 
works 

m = Kf (2.19) 

The direct problem in radar imaging refers to the mapping from the target to the 
quantities that can be measured by the radar. Consequently, the inverse problem is 
concerned with ‘reproducing’ the original target from given data and knowledge of the 
direct problem (in the case of radar imaging, it will be radar measurements and 
knowledge of the scattering model). 

1. Well-Posed and Ill-Posed Problems 

One of the issues for inverse problems is their ill-posed nature. An ill-posed 
problem is one whose solution is not unique and/or does not exist for any data and/or 
does not depend continuously on the data. In the case of a band-limited system, the 
solution of the inverse problem is not unique. This is because the imaging system does 
not transmit information about the target at frequencies outside the band of the measuring 
system. 


20 



Consider discretized measurement systems in which the measurements are 
collected at specific points, and so Equation (2.19) can be written in matrix form. And 
the basic goal of getting a good estimate / from the measurements m can be notionally 
represented as 

/ = K l m ( 2 . 20 ) 

It can be observed that k , which represents the measurement system, imposes certain 
properties on the measurement system. In general, the matrix k is NxM(M>N) as 
the object function is infinitely dimensional and will not be confined to the limitations of 
the measurement system. Consequently, Equation (2.20) has more unknowns than there 
are linearly independent equations and so the system cannot have a unique solution. 

Since the object space (M-dimensional) is typically larger than the measurement 
space (N-dimensional), there exists what is known as the nullspace which consists of all 
vectors / - icf or, physically, all the things that cannot be measured. Nullspace 
accounts for measurement artifacts and is determined by the kernel k . 

The accepted approach for solving inverse problems which are ill-posed is to 
search for approximate solutions satisfying additional constraints based on the physics of 
the problem [22]. This set of approximate solutions corresponding to the same data 
function is the set of objects with images close to the measured one and is expressed in 
what are known as the normal equations 

K T Kf = K T m (2.21) 


2. Data Reconstruction - Regularization 

Equation (2.21) specifies the set of conditions that must be satisfied by the least- 
squares solution 7 which can be rewritten as 

7 = (K r K)- l K T m ( 2 . 22 ) 

Recall that the measured data suffers from finite dimension and noise contamination. 
These are problems since ( k t k ) is poorly conditioned and the effect of noise is 
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magnified (i.e. errors increase) as the dimensionality of the problem gets larger. Two 
techniques will be highlighted here, namely the singular value decomposition (SVD) and 
Tikhonov regularization, to mitigate the effects of noise when estimating / . 

From linear algebra methods, the matrix (k'k) 'k t can be expressed in terms of 
its SVD [19]: 

(k t k)- 1 k t = UDV t (2.23) 

where the matrices U, V contain the eigenvectors of (k t k), and D will have diagonal 
elements of the form 1/ ^A. (where A. are the corresponding, non-zero eigenvalues). A 

consequence of the kernel k being bounded is that the eigenvalues necessarily form a 
sequence that gets arbitrarily close to zero as the number of measurements increase. The 
diagonal elements of D get large when the higher dimensional data are used, and if the 
data are contaminated by noise, from Equations (2.22 and 2.23) that noise will also be 
multiplied by a large value. Hence, one simple approach to mitigate the noise effect is to 
truncate the SVD representation when A. falls below some chosen threshold value. The 

corresponding diagonal element is set to zero and this method is known as the truncation 
filter. 

Tikhonov regularization modifies (i k t kY x k t to a better conditioned form by 
adding a fixed scalar a to the original normal equations. Equation (2.22) can then be 
rewritten in the form 

/ = (k t k + aI)~ l K T m (2.24) 

When a = 0, the original normal equations are returned, so it can be expected that small 
values of a will not change the system description too much. 

After introducing a , the diagonal elements of D in Equation (2.23) become 

D, = (2.25) 

A. + ct 

i 
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When A.? a , D « 1 / JJ~. (as before). When A, = a, however, D « / a and will 

not become arbitrarily large. Consequently, the noise problem is controlled. 
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III. IMAGING DATA MODEL 


This chapter develops an approach to imaging moving objects from scattered 
waves. Many of the ideas and concepts introduced in the preceding chapter will be 
extended to the three-dimensional case in order to develop the imaging model. 

A. INTEGRAL EQUATION APPROACH TO SCATTERING 

We have seen in the preceding chapter that for the one-dimensional scattering 
problem from a moving plate, the signal that returns to the radar is a time-delayed and 
frequency-shifted version of the transmitted signal. This is also the case for three- 
dimensional scattering from a moving, point-like target. A point target is the traditional 
model used for developing imaging analysis and will also be used for the thesis. The 
representation of a target as points is appropriate because for high-frequency waveforms 
scattered from smooth conducting surfaces, the net scattering effects from the surface is 
equivalent to the contribution from the entire series of points known as specular points 

[19]. 

Applying traditional electromagnetic approaches to scattering will typically 
involve solving a differential wave equation and imposition of appropriate boundary 
conditions at the surface of the scatterer. In practice, this process is difficult to 
implement for there is no coordinate system that can represent the boundary conditions 
for complex targets we encounter in practical scattering applications. 

Hence, it is convenient to reformulate the scattering problem. It turns out that 
there exists an equivalent integral equation for every differential equation, and this 
integral equation can intrinsically incorporate any boundary conditions in the solution. 

We take the scattering body to be located in free space and describable by an 
index of refraction n(x) . Consider a monochromatic wave, with harmonic time 
dependence (any signal can be built from these individual frequency component via 
superposition). The reduced wave equation for the field (Helmholtz equation) is 
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V 2x F(x) + k 2 n 2 (x)W(x) = 0 (3.1) 

where k is the wave number in free space (note that n(x) ^ 1 on the scatterer). Splitting 
the field 'P into the incident field ¥ representing the incoming wave, and the 


scattered field x V v . ri „ and simplifying. Equation (3.1) yields 

- nHx))(v,Jx) + V_(*)) 

where the scattering density p(x) = k 2 (l-/7 2 (v)) is a source factor completely 

determined by the index of refraction of the scatterer and is zero outside of the scatterer. 

A solution to Equation (3.2) can be obtained based on the application of Green 
functions , and we can write the solution as 

(3.3) 

This is a Lippmann-Schwinger equation. It can be observed that the problem of solving 
equation (1.1) has been reduced to finding the Green function Gfx'.x) in free space (k 
denotes a specific frequency). The scattered field is also written as an integral over the 
scatterer D since p(x) = 0 otherwise, and hence automatically satisfies the required 
boundary conditions for a known target. 

It is appropriate to make several additional observations with respect to the 
Lippmann-Schwinger equation: 

(1) There are generally many functions that can serve as Green functions 
Gf x'.x) in free space. While there are no bounding surfaces (other than those defined 
by the scattering object), we are interested in G Ax', x) that obey the Sommerfeld 

radiation condition Y\mG k (x' ,x) oc - 
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(2) The inverse scattering problem - i.e. that of determining p given and 
¥, can also be approached using the Lippmann-Schwinger equation as a model for 
'F . 

scatt 


B. LINEARIZED DATA MODEL (TIME-VARYING SYSTEMS) 


Enforcing the Bom (“weak scatterer”) approximation, where the object 
distribution (points) is not self-interacting, multiple scattering events are treated to be 
negligible in comparison with the primary scattering events. The resulting linearized 
(approximate) data model is 

'f',./*) = lllp,wvp<rw,jxv 3 x' (3.4) 


This solution is one for the Helmholtz equation. For time-varying systems, the 
time-dependent wave equation and the time-domain Green function will have to be 
incorporated within the model 


G k (x\x)^> g(x\x;t\t) 


S(t - 1'— \ x'- x\/c) 
4n I x'— x I 


(3.5) 


For a signal s (t ) transmitted from position y at starting time —T , the incident field at x 
is 


¥ inc M = 


s y (t + T y -\x - y\lc) 
4n I x - y I 


(3.6) 


and so the scattered field y/ scatt (y,z,t) at time t and position z obeys (Bom 
approximation) is 


¥ s , 


r(y,Z,0= JJ 


8{t-t'- \z-x'\!c) &(t' + T y -\x'- y l/c) ^ )w3 


4n\z~ x'\ 


4n\x'- y I 


p(x')d x'dt' (3.7) 


where SL denotes the second time derivative of ,v (l) and enters as a consequence of the 
second time derivative in the wave equation. 
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c. 


REFLECTIVITY FUNCTION FOR MOVING TARGETS 


The reflectivity function p(x) or scattering density first introduced in Equation 
(3.2) is completely specified by the target scatterer. It is a scale factor of the received 
signal strength and an important characteristic of radar imaging. For a moving target, a 
time-varying reflectivity function p(x,t ) will be used. For an imaging scene containing 
multiple moving targets (or scatterers), the scattering model can be modified as follows: 
Fet p v {x-vt)d 3 xd 3 v be the scatterers in the volume element d 3 xd 3 v of phase space 
centered at position x and velocity v. Choose T to be such that the scattering density 
p = p (x) describes those targets moving with velocity v at time 1=0 (i.e., the 
transmitting antenna is activated at time t = -T ). The spatial scatterer density centered 
at time t and position x is 

p(x) = ^ p v (x-vt)d 3 v (3.8) 

and the scattered field obeys 




,(y,z,t) = JJJ 


S(t-t'-\z-x'\/c) 

An I z - x' I 
'+ r - I x'— y I Ic) 
An I x'- y I 


(3.9) 


-p v (x'- vt')d 3 vd 3 x'dt' 


Making the change of variables x' —» x = x'— vt' (i.e., change to a frame of reference in 
which the scatterer p is fixed) to obtain 




S(t — t'— \ x + vt'— z\lc ) 

An I x + vt'- z I 

T - I x + vt'- y I Ic) 

x —--- 

An I x + vt'- y I 


p v (x)d 3 vd 3 xdt' 


(3.10) 


The physical interpretation is that the wave that emanates from y at time -T encounters 
a target at time t'. This target, during the interval [0,C], has moved from x to x + vt'. 
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The wave scatters with strength p v {x) and then propagates from position x + vt' to z, 
arriving at time t. 

Simplifying the notation, and writing 

R t T (t) = x + vt - z, R =\ R \, R=R/R 

then 




t (y,z,f) = JJJ 


8{t-t'~ R x z (t')/ c) 
4 ttR XO 


'+T-R (t')/c) 

---—- p (x)d vd xdl ' 

4 nR (t') 

x,y v 7 


(3.11) 


The variable t’ in the argument of the delta function appears in several terms making it 
difficult to carry out the integral over t’. To address this, let t' = t xv (t ) denote the 

implicit solution of 

t-t'~ R xz (t')/ c = 0, i.e., t-T v (t)-R xz (t xv (t))/c = 0 


Then, performing the integration over V yields 

Vscan^^) = \\ 


(t (t) + T -R (;t (t))/c) 

\ x,v v ' y x,y v x,v v / 


p [ (x)d\’d\x 


(4^)XAvW)*,Av«) 


Substituting t x v (t) = t - R x Xt x v (t)) / c for the first occurrence of 1 (t) in this last 
equation allows for a more symmetric version 




# (y,z,0 = JJ 




t + T -R (t ( t))/c-R (t (; t))/c 

y x,z v x,v v 77 x,y v x,v v 77 


) 


< 4 *>' KSK.S»)R .,(*„(») 


p (x)d vd x (3.12) 


D. FURTHER SIMPLIFYING APPROXIMATIONS 

Equation (3.12) applies to very general circumstances. It applies to rapidly 
moving targets (where the start-stop approximation is not valid), and is appropriate for 
arbitrary transmitted waveforms. 
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However, Equation (3.12) is also a manifestation of retarded-time problems, 
which are invariably complex. The complications arise as a consequence of the 
definition of retarded time t x Jt ), which is obtained implicitly from the relation 

t (t) = t-R (t (t)) / c . A more tractable result can be obtained when the scatterers 

X,V x 7 X,z x,v v 77 

are known to be moving in such a way that R Jt) can be expanded in a Taylor series 

around t=0 and approximated by retaining only the terms linear in/- i.e., when the 
scatterer is “slow moving”. 


1. Slow-Mover Approximation 


Assume that (I v If) and (Ivl 2 t 2 x co / c) are much less than I jc — z I and 

v 7 v max 7 

I x - y I (where co nm denotes the maximum angular frequency of the transmitted signal 
5 ). In this case, 

y ’ ’ 


R x Jt) =1 z-(x- vt) 1= R xz ( 0) + R xz (0)Q>t + ... 

where R x JO) = x - z , R x .(0) =1 R _(0) I, and R x ^(0) = R _(0) / R x _(0) • Substituting 
this result into the definition of retarded time yields 

EW«>-(*„(0) + J»„(0)g ’Ijoyc 

l-RJ0)lc 
l + R x J0)Q>/c 


Inserting this approximation into the result for y/ ( y,z,t) yields 


¥ 

’ St 


( (y,z,0 = JJ 


<(a 


I a 

y V 


t-R (0 )/c -R ,(0 )/c + T 


x,y 


1 


pXx)d i vd 3 x (3.13) 


(4/r) R xz (0)R xy (0) 


where 


l-ig ro ,(0)9>/c 
1 + R xz (0)q>/c 


(3.14) 
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is the Doppler scale factor. Similar to the traditional Doppler shift, it is closely related to 
the velocity component in target-transmitter direction. 


2. Slow-Mover and Narrow-Band Approximation 

Typical systems are narrow-band for which the transmitted signal is of the form 


where §/c(t) is slowly varying (as a function of t) in comparison with e “° y (where co is 
the carrier frequency for the transmitter at position y). In this case, the time derivatives of 
s are dominated by the e "° yt factor, and £x(t) « -ci) 2 Sj(t)e my '. Moreover, since §/<(t) 
is slowly varying, further approximations can be made 


,?/c( a x t)e 


(3.15) 



3. Slow-Mover, Narrow-Band and Far-Field Approximation 

When the transmitter-to-target and target-to-receiver distances are large in 
comparison with the scene dimensions, then \x + vt\ and I x + vt\ 2 xoo /c can be 

L max 

assumed to be much less than either I z I or I y I. And the expansion 

R x z {t) =1 z-(x- vt) 1=1 z I -ixjv + vt) +... 
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can be applied (and similarly for R v (f)). Substituting the expansions into previous 
approximations (Equation (3.16)) then yields 


¥ sca ly^ t) 


a y \\e 9 *'e~ ia,ya * v> 

(4tt) 2 I - II y I 

x + T y - (I z I -itpr+1 y I — tagr) / c y (x)d 3 vd 3 x 


(3.18) 


where now 


cv =co 

' x,v y 


I y I -pgr - cT y + a v (I z I -PCpr) 


/ c 


(3.19) 


Note that in Equation (3.16), the expansions for R (t), R (?) are used differently. In 
the amplitude factor the approximation R _(?) =1 z - (x + vt) l»l z I holds since 

I z I? I fig.i' + vt) I; in the phase factor, however, the product (pgr x a> max / c) can still be a 

large fraction of 2 n , and must be retained. In addition, the signal envelope §/ can still 
vary significantly over I z I -pgr despite being slowly varying, hence the first order term 
is retained. 

The Doppler scale factor also reduces to 

a v *« 1 ■- (3»+ m / c (3.20) 

u 1+pCp/c 

where the binomial expansion is used since lvl/c= 1. The quantity 0 ) x /) is the 

Doppler shift (where /? = -( y>+ p)g» / c) and is observed to be dependent on the 
“bistatic” vector ( y>+ p). 

Inserting T =1 y I /c and setting k = co / c 

cp x v =k y lzl-^ y (p+p)^x + (pgz-x))v/c] (3.21) 


32 



E. IMAGING VIA A FILTERED ADJOINT 


In developing the correlation receiver, it was observed that the “fit” between two 
(complex-valued) functions is given by the cross correlation. For correlation reception, 
the correlation integral seeks out the component of received signal that “matches” the 
time-delayed, frequency-shifted version of the transmitted signal. Large values of the 
correlation integral indicate a strong resemblance, while small values indicate weak 
resemblance. The values of the time-delay and Doppler shift parameters, which 
maximized the correlation integral, were also the parameters most likely to represent the 
target. 

Using the same notion, imaging can be achieved by determining the position and 
velocity parameters that represent the best “fit” to the measured data. These will also be 
the parameters that localize the position and velocity of the unknown scatterers in phase 
space (position-velocity) space. Cross correlation in the time domain, is multiplication in 
the Fourier transform domain ( Fourier correlation theorem) - hence the association with 
filtering since filters are multiplication operators in the Fourier transform domain. 

The model developed for i// (y,z,t) depends on the position of the scatterer 

jkt = (jCj,jc 2 ,jc 3 ) and its velocity v = (v 15 v 2 ,v 3 ) . Using this form for y/ scatt (y,z,t) to build a 
six-parameter scattering model based on the arbitrary position p = (p r p 2 , p 3 ) and 
arbitrary velocity u = (u^ip.u,) , the objective will be to determine the values of p and u 
that maximize the cross correlation between the arbitrary model and the measured data. 

In terms of the Doppler shift, the scattered field (from known scatterers p v (x) 
with positions x and velocities v) is 

2 -io (t-\z\/c) 

—CO 6 m /• r A 

Vscanty^ ~ ~\ 2 , „ , 11 ex P M-X - v(f - zc(z - X) / c)\ } 

(4;r) 2 1ZII y I 1 y } 

x §/o(t - (I z I -(jj-T p)C|r) / c\p v (x)d 3 vd 3 x 
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The effect of data shifted in time is an image that is translated in range. Since this 
has no bearing on the imaging process, the scattered field result can be simplified by 
substituting t' = t-\z\/c which yields 

-are " 0y> rr r 

Yscanty^') * * | ^ H | JJ { ~ik y (y>+ ~ V(t'+ HOC / c)] ] 

x §Zo(t'+ (y>+ j3)gr / c\p v {x)d\d i x 

Ignoring the intensity prefactor, a parametric model representing the field from an 
unknown point source is 

f y(y,z,t ') = ,<v exp i-ik ( y>+ - u(t'+ sq p / c)]) 

/ y J (3.22) 

x §Ah'+(y>+ v)Q/) / c) 

And the image is created as follows: For each p and u, determine the function 

I(p,u) = JJJ i/f scatt (y,z,t')V'*(y,z,t')dt'd m yd n z (3.23) 

where m, n depends on the configuration of the transmitter(s)/receiver(s), and the 
integrals are over all values of t ', y and z for which measured data is available. 

I(p,u) = JJJ i// scatt (y,z,t')i//\y,z,t')Q(o),t',p,u,y,z)dt'd m yd n z 

t / x fff/o/ \ k °r * v (o+e)J/>-«(r+sig/7/f)] 

/(/?,«) =-JJJ Q(co,t ,p,u,y,z)e e y 

x .?/o (t'+(p+ v)Qp / c(y,z.,t')dt'd m yd"z. 
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IV. IMAGING POINT SPREAD FUNCTION 


This chapter begins by looking at the radar ambiguity function, and the 
relationship between the transmitted waveform, and the accuracy with which target range 
and velocity can be estimated. The imaging point-spread function (PSF) is also derived 
based on the imaging scheme developed in Chapter III. 

A. RADAR AMBIGUITY FUNCTION 


For a stationary target, the output of the correlation receiver is the cross 
correlation between (1) the received signal plus noise and (2) the transmitted signal. In 
many radar applications, however, the target is moving so that its echo signal has a 
Doppler frequency shift. The output is hence a cross correlation between the received 
signal and the Doppler-shifted transmitted signal. 

The nature of the correlation receiver output as a function of both time and 
Doppler frequency is important for understanding the properties of a radar waveform, in 
particular its effect on measurement accuracy, target resolution, and ambiguities in range 
and radial velocity. These aspects of the correlation receiver output will be examined in 
turn. 


When the received echo signal is large compared to noise, the output of the 
correlation receiver may be written as the following (see Equations (2.15) and (2.16)) 


7 (y,r) = JJ p(u',T')x(v-u\T-T')e 


i-(u+v')(T-z') 

2 dr'do' 


where %(u - v ', r - r') is the radar ambiguity function. The problem of radar imaging is 
concerned with estimating the location and strength of the point scatterers that are 
assumed to make up the target. Consequently, understanding % is essential as the end 
state involves ‘inverting’ the radar data model. By definition, we have 


Z(o-o',r- r') 


1 







gi(v-u')t' 


dt' 


(4.1) 


—co 
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1. 


Basic Properties 


For simplicity, the origin is chosen to be the true target time-delay and frequency- 
shift. This assumption will be carried throughout the rest of the discussion. Some key 
properties that follow from Equation (4.1) [Refs. 14, 18 and 19] are listed 

Signal energy : \%(v, r)| < | j(0,0)| (4.2) 

oo 

Ambiguity volume : JJ |j(t>,r)| dzdv = 1 (normalized signal) (4.3) 

—oo 

Symmetry : \%(-v,-z)\ = I %(o, r)| (4.4) 


Frequency - domain: %(o , r) 



— v)s* (co + — o)e ,COT dco 
2 mc 2 


(4.5) 


Equation (4.2) states that the maximum value of the ambiguity function occurs at the 
origin (true target time-delay and frequency-shift). Equation (4.3) asserts that the total 
volume under the ambiguity surface is a constant. This gives rise to a radar uncertainty 
principle: choosing a signal s(t) so that the ambiguity surface will be narrow in one 
dimension will cause it to be correspondingly wide in the other dimension. In the case of 
range-Doppler imaging, a waveform with good range resolution has poor Doppler 
resolution and vice versa. 


2. Range Resolution 


Equation (4.5) presents an intuitive look at range resolution. Consider a fixed 
target (o = 0) and rewrite Equation (4.5) as 


lj(0,r) 1= 

2 n 


[is. (co) \ 2 e^ T doj 

j inc v ' 


(4.6) 


Recognize the form of Equation (4.6) as the inverse Fourier transform (of the transmitted 
power spectral density). Then, in order for I j(0, r) I to have high resolution in z 
measurements (ideally a delta function), the power spectral density has to be identically 
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one; such a signal will have infinite energy or broadly supported in the Fourier domain. 
Hence, it is clear that better range resolution is obtained when the signal bandwidth is 
broad. 

A more rigorous study for estimation errors (in the presence of Gaussian noise) is 
developed in [1] by examining the behavior of the ambiguity function in the 
neighborhood of its main peak. This study shows that estimation error of v and z is 
described by an ellipse, and that time-domain resolution is inversely proportional to 
frequency-domain bandwidth. Hence, bandwidth is inherently important to problems in 
radar target imaging. A further interesting result makes use of the definition of 
bandwidth in terms of the second moments (probabilistic variance) of the power spectrum 
to yield (for a narrow-band signal s(t) = a(t)e mn ). We can write 



da(t') 

dt' 


dt' + 


1 U d®(t') Y 
2II dt' ) 


I ait') I 2 dt'-a>l 


(4.7) 


where (3 represents bandwidth. The useful relationship illustrates a nonlinear O(r) 
(phase modulation) will increase the bandwidth in comparison with signals that are 
modulated in amplitude only. This is the basis for pulse compression (see Chapter V), a 
technique of practical importance since it allows for the creation of fine range-resolution 
waveforms that are also of long duration (with sufficient energy on target). 


3. Doppler Resolution 


To study Doppler resolution, consider a target whose range is known. Thus, 
setting r = 0 in Equation (4.1) yields 


\%(o,0) 1 = 


[ I s. (O I 

J inc v ' 


2 ivt 


dt' 


(4.8) 


It is evident that better Doppler resolution is obtained from a long duration signal. This is 
in stark contrast to range resolution requirements. In addition, the Doppler resolution is 
determined only by the amplitude modulation of the signal and not by the phase. 
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B. IMAGE ANALYSIS 


From Chapter III, an explicit form for the image is 

I(p,u) = -JJJ 

X S/c(t'+(p+ ©)gp / c)y/ scan (y,z,t')dt'd"'yd"z 


(4.9) 


To examine the performance of the imaging scheme, the ideal (expected) data are 
inserted back into the image equation. Substituting for y/ , (y.z.l') yields 


'<*«) = M 


co 2 y Q(co,t\p,u,y,z ) 
(4;r) 2 1 z II y I 


pM) 


exp Jz7c y (jj+ i))(Jp - x-(u- v)t'~ u(sQj) / c + v(ig*:) / c] J 
Wo(t '+ ( jj+ xj)Q/) / c^/d(t'+ (p+ £i)C(r / c')cl\’d i xdt'd" l yd n z 


(4.10) 


The distances I z I and I y I arc scaling amplitude factors and not essential for 
imaging purposes. Thus, choose 


ru - \ ( 4 ^) 2 I Z II y I r, , 

Q{co,t ,p,u,y,z) = -j—— J(p,u,y,z) 


(4.11) 


CO 


where J(p,u,y,z ) depends on the geometry and is chosen to compensate for a Jacobian 
that results from the integral of the variable t’ 

I(p,u) = JJJ JJ$p(f'+ ( j»+ ©)gp / c)?/o(T+ ( y>+ ©)gr / c)e~^ (y+z)! *““ v) ' 

x exp ^k y (p+ V)Qp - x - u(zQp) I c + v(cgc) / c] J- (4.12) 

x J(p,u,y , z)p v (x)dt'd 3 vd 3 xd m yd"z 

which is 

I(p,u) = JJ K(p,u,y,z)p v (x)d\d 3 x (4.13) 

where 
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K(p,u,y,z ) = JJ exp ^-ik y (p+ $)Qu(gQp) - v(agv:)] / c j 
x | §/o(t '+ ( 59 + £)gp / c^?/o^'+ ( 59+ cDgv / c)e 
x J(p,u,y,z)d'"yd n z 


-ik ($+i>)c(w-v)r 


dt' 


(4.14) 


is the point-spread function describing the behavior of the imaging system. Apply the 
change of variables t = t'+ “(?+ £>) 2 [p + x) / c and set r = ( 59 + ©) 0 |: p-x) / c to obtain 


K(p,u, y,z ) = JJ exp { -ik y (p+ g)$u(gQp) - v(igr)] / c j 

x exp |i^( 59 + v)Q(p -x) + ^(u- v)( 59 + 0)0 ip + x)J / cj 


(4.15) 


f n *( 1 V/ 1^ -ik (H+v)Qi<-v)t , 

J §/o^t + - rj sjcjj --Tje y dt 


J(p,u,y,z)d'"yd n z 


The integral in square brackets is the radar ambiguity function 


» f 1 A f 1 A 

Z(o, r) = J S^t + - rj S/d^t - - rj 


X is determined by the transmitted signal and I x I attains its maximum value at p = x , 
11 = v (for general configurations of transmitter/receiver). So the imaging PSF can be 
written as 

K(p,u, y,z ) = JJ exp {-/£( 59+ p)Qu(HQ)) - v( ®C|r)] / c j 

x exp jik v ( 59+ &)(J( p -x) + ^(m - v)(j9+ i))G( p + x)J / c| (4.16) 

x x(k y (y>+ ©)S( M - v),(j 9 + ii)g(p - x) / (p,u,y,z)d'”yd"z 
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V. PULSE COMPRESSION 


Detection of signals in noise and resolving capability are both fundamental to 
radar imaging. As discussed in Chapter II and IV, the use of correlation reception and a 
transmitted waveform with non-linear phase modulation provides means to enhance 
performance in detection and resolution. Pulse compression is a technique that fuses 
these concepts and involves transmitting long coded pulses, together with signal 
processing methods such as correlation reception. This is of practical importance since it 
allows for the creation of fine range-resolution waveforms that are also of long duration 
(this increases radiation energy directed on target without relying solely on boosting radar 
transmitted power). This chapter gives a physical interpretation of pulse compression 
and looks at a specific phase-coded signal known as the chirp. 

A. REVISITING CORRELATION RECEPTION 

Recall in Chapter II, the output of the correlation receiver is a correlation integral 
between the received signal and a time-delayed, Doppler-shifted version of the 
transmitted signal. It is illustrative to explain the phenomenon of pulse compression with 
a graphical example. Using the simple one-dimensional imaging case in Equation (2.12), 
the correlation receiver output can be written as 

oo 

77 ( 0 = \s*(t')s rec (t'+t)dt' (5.4) 

-00 

which is a correlation between s and s . If s = s (reflectivity function is 1), Equation 

(5.4) is called an autocorrelation. In Figure 5. the signal at the top represents the 
transmitted waveform. The lower graphs represent the received signal for different shifts, 
according to Equation (5.4). When a signal shifted by t is multiplied by the waveform at 
the top and the product is integrated over t’, the resulting number is plotted at position 1 
on the graph to the right. 


41 





Figure 5. Phenomenon of pulse compression where the energy is concentrated at a 

single delay time [18]. 

Figure 5 shows that, although the output signal has a longer time duration than the 
original signal s, its energy is more concentrated at a single delay time. It is a 
consequence of the Cauchy-Schwarz inequality that the highest peak occurs when the 
signals are not shifted relative to each other. The smaller peaks, called sidelobes, are 
undesirable and an open problem is to design pulses whose autocorrelations have the 
lowest and fewest sidelobes. An issue that has direct bearing on signal detection as the 
received signal is expected to “look” like a time-delayed, Doppler-shifted version of 
itself. 

B. PHASE CODING - CHIRPS 

A radar system makes the most efficient use of power when it transmits a 
constant-amplitude waveform. Consequently, most radar waveforms are designed with 
variations only in the phase of the signal. The most important and most commonly-used 
radar waveform is the chirp. A chirp is a constant-amplitude signal whose instantaneous 
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frequency varies linearly with time; chirps are therefore also called Linearly Frequency 
Modulated (LFM) waveforms. Linear variation of the instantaneous frequency implies 
that dtp / dt = a> mm + yt . The coefficient y is called the (angular) chirp rate. A chirp is 
thus of the form 


s(t) = e 


i(co .t+—yt ) 


M [o ,r]W 


(5.5) 


where m (f) is 1 in the interval [0, 7] and 0 otherwise. 

A chirp with positive chirp slope is called an upchirp; one with a negative chirp 
slope is a downchirp. The instantaneous (angular) frequency of the chirp varies from 
co . to co = co . + yt, with a center frequency of <x>,=co + yt / 2 . Thus the 

instantaneous frequency of the chirp varies over the interval [<y () - Ty / 2,a> 0 + Ty / 2] . 
The power spectrum of a chirp is roughly constant over the frequency band 
[® 0 - Ty / 2,co 0 + Ty / 2], thus the angular frequency bandwidth Ty can be estimated 
by looking at the range of instantaneous frequencies (see Figure 6). 



uo - T-7/2 ljq \ rV/3 


Figure 6. Spectrum for a chirp pulse of length r [1]. 


Chirps are very commonly used in radar imaging. Pulse compression for chirps 
has the following intuitive interpretation: When the pulse is transmitted, different parts of 
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the pulse are coded with different instantaneous frequencies. In the pulse compression 
process, the different frequencies are delayed by different amounts of time, so that all the 
energy emerges at the same time [18]. 

C. AMBIGUITY FUNCTION FOR COMMON WAVEFORMS 


From the preceding discussion, the “ideal” ambiguity function will typically 
consist of a single peak of infinitesimal thickness at the origin and be zero everywhere 
else. It will closely approximate an impulse function, and have no ambiguities in range 
or Doppler frequency. The infinitesimal thickness at the origin will permit the time delay 
and/or frequency to be determined simultaneously to as high a degree of accuracy as 
desired. Such a highly desirable ambiguity function, however, is not theoretically 
allowed. The properties of the ambiguity function implicitly mandate that trade-offs exist 
between the measured dimensions. 

In practice, other than the trade-offs in dimension resolution, waveforms also 
generally have significant response outside the narrow region in the near vicinity of the 
origin, i.e. there are many possible ambiguous (blind) range and Doppler. Ambiguities 
are a consequence of discontinuous waveforms, such as in a pulse train. The resolution 
for single rectangular pulse, single chirp and coherent pulse trains will be discussed next. 


1. Resolution for a Single Rectangular Pulse 


Time domain impulses have infinite bandwidth, hence they are idealized signals. 
More realistic waveforms can be modeled, for example, as 


Jl if 0 < t < T, 

[0 otherwise. 


(5.6) 


for which s. nc (t - 1 / 2z)s* ne (t + 1/2r) is non-zero in the interval t e(l r I /2, T— I r I /2). 


44 



The ambiguity function is then 


00 i i 

r 1 1 r T-lrl/2 . , 

X(v,t)= fs. (f'--T)s* (t'+-T)e w 'dt'= f e wt dt 

7 J m c v 9 ' 9 Jlrl/2 


= —e 
io 


iu(TI 2) 


giv(T /2—lrl/2) _ ^~io(772-lrl/2) 


(5.7) 


= (T- I r l)e" ;r/2 sinc 


1 


y(r- I z-1) 


For u = 0, then the delay resolution (by the peak-to-first-null measure) as 
observed from the ambiguity surface l^(0,r)l 2 , is A t = T (zero-crossing at t = T), 
which corresponds to a range resolution of A R = cT / 2. For r = 0, the ambiguity 
surface I %(o,0) I 2 yields A o = 1 IT since the first zero of the sine function occurs when 
its argument is n. The Doppler shift is related to down-range relative velocity by 
v n = -2vl> 0 / c (see Equation (2.4)), and the corresponding velocity resolution is 

Av = A () / 2 T . 


2. Resolution for a Single Chirp 


This is another example of a ridge type ambiguity function produced by linearly 
frequency modulating a rectangular pulse over a bandwidth ft. The pulse width T is 
large compared to 1 / /?. The frequency modulation increases the spectral bandwidth of 
the pulse so that /3T ? 1. Since the pulse width T and the bandwidth [3 can be chosen 
independent of one another, the time-delay and frequency accuracies are independent of 
the other. 


The chirp signal is one whose frequency changes linearly with time 


s. (t) = < 

inc x ' 


r ,i 

ex p^-*-n~j 


0 


if 0 < t < T, 
otherwise. 


(5.8) 


45 



The constant y is known as the chirp rate. The ambiguity function is then 


co , 

r 1 1 rT- lrl/2 .. 

Z(u, z) = \s (t'~ - z)s* (t'+ - r)e wr dt' = f e ,(u+n) ' 

' J O O J lrl/2 




1 


i(v+n)TI2 


i(o + yz) 

= (T- I r l)e' ( " +rr)r/2 sinc 


^/(D+7r)(r/2-lrl/2) _^-*(D+7T)(772-lrl/2) 


(5.9) 


-(u + ^r)(r-l rl) 


Observe that a chirp leads to an ambiguity function that is a .vine-function in both the time 
and frequency domain. It is instructive to examine this property again by looking at the 
signal resolution. 

For z = 0, similar result of Ao = 1 / T is obtained as in the case of the rectangular 
pulse. This is because the amplitude modulation is identical in both cases. 
Consequently, Doppler resolution is the same. However, with phase modulation, it is 
expected that the delay resolution (hence range resolution) will be enhanced. 

For o = 0, when z > 0 the first null occurs when n = (T - z)yz / 2. Reference 
[19] shows that the expression /z = (T-z)yz/2 has a quadratic equation form with a 
solution Az = 2zz / yT, which corresponds to a range resolution of A R = c/ J3 where 
/? = 2 7zyT is the bandwidth (in Hertz). For r = 0, Ao =1 %(o, 0) l 2 = 1 IT since the first 
zero of the sine function occurs when its argument is tz . The Doppler shift is related to 
down-range relative velocity by o D = -2vo {) / c (see Equation (2.4)), and the 

corresponding velocity resolution Av = A / 2 T . 

The pulse compression ratio is defined as the product of pulse spectral bandwidth 
(] and the uncompressed pulsewidth T. Comparing the range resolution between the 
single rectangular pulse and single chirp, 


A R 


rectangular 


A R 


chirp 


cT ! 2 pT 
c!p ~ 2 


(5.10) 
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In other words, phase modulation improves the range resolution by a factor of one-half 
the time-bandwidth product. 

3. Resolution for Coherent Pulse Trains 


Following the analysis given in [18], a train of N identical pulses is considered: 


s(t) = 


1 

Viv 


N-l 


n =0 


(5.11) 


where T is the pulse repetition interval, which might be on the order of 1ms. The 
reciprocal 1 / T R is called the pulse repetition frequency (PRF). 


The ambiguity function for the pulse train is 


1 N-IiV-l 

ZO>.r) = ±ZZl 

1V n =0 m=0 


u(t + T'- mT )u(t'~ nT )e iUT dz' 


1 

N 


N-l N-l 

nnxX T ~( n - m)T « - °y" mTR 


n =0 m=0 


(5.12) 


where x denotes the ambiguity function for u and where we have used the substitution 
t" = z'~ nT R . A lengthy calculation involving rearrangement of the order of summation 
and the summing of a geometric series results in 


i N-l 

X(u,z) = — X Zj,*-pT R ,D)e uirR(N ~ l ~ m 

™ p=-(N-l) 


sin tzvT r {N- I p I) 
sin( 7 zuT ) 


(5.13) 


If the pulses u are sufficiently well-separated so their ambiguity functions x„ do not 
overlap, then we take the absolute value of Equation (5.13) to obtain 


1 N-l 

\x(v,z)\=— Yj pT R ,o)\ 


N 1 

p=-(N- 1) 


sin 7zvT r (N- I p I) 


sin {7wT ) 


(5.14) 
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Such an ambiguity function has a “bed of nails” appearance (see Figure 7), with peaks at 
r=pT R , p = -(jV-l),-(iV-2),...,0,l,...,(iV-l) and o such that 

tvuT r N = (m +1 / 2) tv , m = 0,1,2,... 


Ix(r,i/)| 2 



Figure 7. Ambiguity function for a train of pulses. 


Following previous analysis, consider the v = 0 cut of the ambiguity surface to 
estimate the range resolution yields 

12f(0,r) 1= -j- 2 I^r-P^KA-Ipl) (5.15) 


where we have used the fact that lim[sinfot]/ sin x-b. We see that the range resolution 

x —> 0 

of the main peak is the same as the range resolution of the individual pulse u, but that we 
now have range ambiguities due to the extraneous peaks at r = pT R . 


To estimate the Doppler resolution, consider the r = 0 cut of Equation (5.14) and 
for well-separated pulses, only the p=0 term contributes 


Yj \Z u (0,o)\ 

A p =-(N- 1) 


si n [ vT r N\ 
sinful 


(5.15) 
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Consider the numerator argument since it varies much faster. It shows that the Doppler 
resolution is 1 / NT , and hence, the Doppler resolution of the pulse train is higher than 
that of a single pulse. 

The many ambiguities produced by a pulse train may appear to lead to a poor 
radar waveform. The time-delay measurement accuracy (determined by pulsewidth T) 
and the frequency accuracy (determined by the pulse repetition interval T ) can be 

selected independently. If pulse repetition interval T R is such that no radar echoes are 
expected with a time delay greater than T R , and no Doppler-frequency shifts are expected 
greater than 1 /T , then the effective ambiguity function reduces to just a single spike at 
the origin whose dimensions are determined by T and T R . In practice, ambiguities can 
also be resolved with different PRF, and many radars employ this type of waveform. 
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VI. RESULTS AND CONCLUSION 


Here, analysis of the imaging point-spread function (PSF) (Equation (4.16)) is 
carried out using Matlab. Equation (4.16) does not restrict the choice of the transmitted 
waveform nor the transmitter/receiver configuration. The goal here is to study the 
behavior of the PSF under different waveforms and transmitter/receiver configuration. 

Figure 8 shows the transmitter/receiver geometries for the simulation performed. 
The receivers are denoted by “zl” to “z8” while the transmitter is denoted by “yl”. For 
configurations of two to four receivers, the geometry as shown in Figure 8(a) is used, 
while Figure 8(b) shows the geometry for eight receivers. 



x-component 

(position) 


(a) 
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x-component 

(posilion) 


(b) 

Figure 8. Transmitter/Receiver configuration for 2, 3, 4 and 8 receivers. 

A. IMAGING PSF FOR SINGLE RECTANGULAR PULSE 

The transmitted waveform is single rectangular pulse defined in Equation (5.6). 
Its associated ambiguity function (Equation (5.7)) is substituted into the imaging PSF, 
Equation (4.16) and implemented in Matlab. Equation (4.16) is capable of determining a 
three-dimensional distribution of positions in space and the corresponding three- 
dimensional velocities (i.e. a six-dimensional phase space image). However, to study the 
effects of transmitter/receiver geometry and translation in position/velocity on the 
imaging PSF, separate two-dimensional plots of the PSF in position and velocity space 
will suffice. 
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1 . 


Localization of Point-Scatterer in Position Space 


The point scatterer is located at position (x v x 2 ) = (0,0). Two-dimensional plots 

of the imaging PSF in position space for two, three, four and eight receivers are given in 
Figures 9(a)-(d) respectively. Clearly, the point scatterers are localized in position space 
about the scatterer position (0,0). However, under the different geometries, the point 
scatterer is represented differently, with the ‘two receivers’ case resembling the least like 
a point. In Figure 9(a), the point scatterer appears longish and rectangular in shape. 
These distortions of the point scatterer are responsible for the image artifacts encountered 
during imaging. The extent of the PSF distortions will determine the level of artifacts 
produced. In Figure 9(d), the PSF appears as a blob and one will expect to see smudges 
in the resulting image. It can be observed that the geometry of the transmitters/receiver 
plays a role in shaping the PSF. 
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Figure 9. 2-D Localization of point-scatterer in position space (for rectangular pulse). 

Figure 10 shows the point-scatterer is relocated to position (x v x 2 ) = (5,5). The 

resulting imaging PSF becomes localized now at the new position (5,5). The PSF retains 
its shape, and the characteristic of translation (in position) invariance for the PSF is 
necessary for a good imaging system. 





x-component 

(velocity) 


Figure 10. Position invariance of imaging PSF. 


2. Localization of Point-Scatterer in Velocity Space 

The point scatterer is located at velocity (v p v 2 ) = (0,0). Two-dimensional plots 
of the imaging PSF in velocity space for two, three, four and eight receivers are given in 
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Figures ll(a)-(d) respectively. Clearly, the point scatterers are also localized in velocity 
space, and similar to the case in position space, the geometry of the transmitters/receiver 
plays a role in the shape of the PSF. The velocity resolution is visibly poorer. 



Figure 11. 2-D Localization of point-scatterer in velocity space (for rectangular pulse). 

When the point-scatterer is moving with a new velocity of (v t ,v 2 ) = (-35,-50), 

Figure 12 shows that the resulting imaging PSF becomes localized now at the new 
velocity (-35,-50). The imaging PSF demonstrates similar translation invariance in the 
velocity space. 
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(velocity) 


Figure 12. Velocity invariance of imaging PSF. 

B. IMAGING PSF FOR SINGLE CHIRP 

In Chapter V, it was shown that a chirp signal improves range resolution. The 
associated ambiguity function for a single chirp (Equation (5.9)) is substituted into the 
imaging PSF, Equation (4.16) and implemented in Matlab. 

1. Localization of Point-Scatterer in Position Space 

Again, the point scatterer is located at position (x t ,jc 0 ) = (0,0). Two-dimensional 

plots of the imaging PSF in position space for two, three, four and eight receivers are 
given in Figures 13(a)-(d) respectively. The PSF is visibly better localized as compared 
to transmitting a rectangular pulse. 
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C. CONCLUSION 

Image artifacts result when motion in the imaging scene is not addressed. In 
recent years, a number of attempts to develop imaging techniques that can handle moving 
objects have been proposed, however, all of these techniques rely on making the 
approximation that a target in motion is assumed to be momentarily stationary while it is 
being interrogated by a radar pulse. This thesis is a study of the new imaging approach 
developed by Cheney, Borden [2] that can accommodate target motion during the 
imaging process. 

In this thesis, the subject of radar imaging from scattered waves is explored and 
applied to derive the corresponding imaging point-spread function for the new imaging 
approach. The simulation result (using Matlab) obtained showed that the imaging point- 
spread function is well behaved, specifically it localizes the target in phase space and is 
translation invariant. In addition, it is apparent that the geometry of the transmitters and 
receivers disposition affects the behavior of the point-spread function, and hence the 
imaging system. 

This work represents the early stages of development of a full imaging algorithm 
applicable to moving targets in a multistatic environment. Follow-on work will include 

further simulation efforts. This thesis used single pulse radar waveforms and a single 
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transmitter, but in reality, fielded radar systems typically utilize pulse compression 
techniques and send out long, coded pulse trains. Neither need the multistatic radar 
system be restricted to one transmitter. In this case, the imaging point-spread function 
has to be adjusted to incorporate more realistic radar waveforms and optimize the 
transmitter/receiver geometries for implementing the imaging algorithm. These will be 
used to support the development of the eventual imaging algorithm. Finally, real world 
target data can then be applied to the developed imaging algorithm to assess its 
performance. 
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APPENDIX: MATLAB CODES 


0,Q,O.Q,O.Q,O.Q,0,Q,O.Q,O.Q,O.Q,O. 

•O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O 


% This program evaluates the Point Spread Function (PSF) 

% The program uses the approach, techniques on imaging from 
Borden, % Brett 


Professor 


% It is a six-parameter scattering model based on arbitrary 
% p=(p_x,p_y,p_z) & velocity u=(u_x,u_y,u_z) 


position 


clear all 


%% Variables to be specified 


x=(10 A 3)*[0;0 ; 0] ; 

% 

actual target 

v=[0;0;0] ; 

% 

actual target velocity 

yl=[2;0;0]; 

% 

transmitter_l 

zl=[-5;5;0] ; 

% 

receiver_l 

z2=[15;15; 0] ; 

% 

receiver_2 

z3=[-2;-15;0] ; 

% 

receiver_3 


z4=[-10;10;0]; 
z5=[0;14.1;0] ; 
z 6 = [ 10 ; 10 ; 0 ] ; 
z 7 = [ 14.1; 0 ; 0 ] ; 
z8=[10;-10;0]; 
u=[0 ; 0 ; 0] ; 
pulse=10 A (-4); 


%% Parameters Defined 

Yl=(1/sqrt(yl(1) A 2+yl(2) A 2+yl(3) A 2)),*yl; 
Zl=(l/sqrt(zl(l) A 2+zl(2) A 2+zl(3) A 2) ) ,*zl; 
Z2=(l/sqrt(z2(l) A 2+z2(2) A 2+z2(3) A 2) ) . *z2; 

Z3= (1/sqrt (z3 (1) A 2+z3 (2)"2 + z3 (3) / '2)).*z3; 

Z4=(1/sqrt(z4(1)^2+z4(2)^2+z4(3)^2)),*z4; 
Z5=(l/sqrt(z5(l) A 2+z5(2) A 2+z5(3) A 2)).*z5; 
Z6=(l/sqrt(z6(l) A 2+z6(2) A 2+z6(3) A 2)).*z6; 

Z7= (1/sqrt (z7(l) / '2 + z7 (2 )^2 + z7 (3)^2) ) ,*z7; 

Z8 = (1/sqrt (z8 (1) / '2 + z8 (2) ^2 + z8 (3) ^2) ) ,*z8; 

J=l; % Jacobian 

c=3 *10^8; 

k_y=2*pi*(10 A 10)/c; % wave vector @f=10 GHz 


%% PSF Range Defined 
p_x_start=-2 0 ; 
p_x_end=2 0 ; 
p_y_start=-2 0 ; 
p_y_end=2 0 ; 

res_p=0.1; % resolution of p 

p_x=p_x_start:res_p:p_x_end; % define vector of x-axis values 
p_y=p_y_start:res_p:p_y_end; % define vector of y-axis values 

%% PSF Velocity Defined 
u_x_start=-200; 
u_x_end=2 00; 
u_y_start=-150; 
u_y_end=150; 
res_u=l; 

u_x=u_x_start:res_u:u_x_end; 
u_y=u_y_start:res_u:u_y_end; 
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Localization in Position-Space 



nx=length(p_x); 
ny=length(p_y); 
K=zeros(ny,nx); 


% length of x = columns in K 
% length of y = rows in K 
% initialize K matrix 


Kl=zeros(ny,nx); 

K2=zeros(ny,nx); 

K3=zeros(nx,ny); 

K4=zeros(nx,ny); 

K5=zeros(nx,ny); 

K6=zeros(nx,ny); 

K7=zeros(nx,ny); 

K8=zeros(nx,ny); 

for m=l:ny 

for n=l:nx 

p=(10^3)*[p_x_start+((n-1)*res_p);p_y_end-((m-1)*res_p);0] 


% for receiver_l 

argl_l=((Yl+Zl). 1 )*((u*((Z1. 1 )*p))-(v*((Z1.')*x))); 
phasel_l=exp(-i*k_y*argl_l/c); 

arg2_l=((Yl+Zl). 1 )*((p-x)+(0.5*(u-v)*(((Yl+Zl). 1 )*(p+x)))) 
phase2_l=exp(i*k_y*arg2_l/c); 
tau_l=l/c*(((Yl+Zl). 1 )*(p-x)); 
neu_l=k_y*(((Yl+Zl). 1 )*(u-v)); 

%rectangular pulse 

amb_l=pi*(pulse-abs(tau_l))*exp(- 


i*0.5*neu_l*pulse)*sinc((0.5/pi)*(-neu_l)*(pulse-abs(tau_l))); 
K1(m,n)=phasel_l*phase2_l*amb_l*J; 


% for receiver_2 

argl_2=((Y1+Z2). 1 )*((u*((Z2. 1 )*p))-(v*((Z2.')*x))); 
phasel_2=exp(-i*k_y*argl_2/c) ; 

arg2_2=((Yl+Zl). 1 )*((p-x)+(0.5*(u-v)*(((Y1+Z2).')*(p+x)))) 

phase2_2=exp(i*k_y*arg2_2/c); 

tau_2=l/c*(((Y1+Z2). 1 )*(p-x)); 

neu_2=k_y*(((Y1+Z2). 1 )*(u-v)); 

amb_2=pi*(pulse-abs(tau_2))*exp(- 


i*0.5*neu_2 *pulse)*sinc((0.5/pi)*(-neu_2)*(pulse-abs(tau_2))); 
K2(m,n)=phasel_2 *phase2_2 *amb_2 *J; 


% for receiver_3 

argl_3=((Y1+Z3). 1 )*((u*((Z3. 1 )*p))-(v*((Z3. 1 )*x))); 
phasel_3=exp(-i*k_y*argl_3/c); 

arg2_3=((Y1+Z3). 1 )*((p-x)+(0.5*(u-v)*(((Y1+Z3).')*(p+x)))) 

phase2_3=exp(i*k_y*arg2_3/c); 

tau_3=l/c*(((Y1+Z3). 1 )*(p-x)); 

neu_3=k_y*(((Y1+Z3). 1 )*(u-v)); 

amb_3=pi*(pulse-abs(tau_3))*exp(- 


i*0.5*neu_3 *pulse)*sinc((0.5/pi)*(-neu_3)*(pulse-abs(tau_3))); 


K3(m,n)=phasel_3 *phase2_3 *amb_3 *J; 

% for receiver_4 

argl_4=((Y1+Z4).')*((u*((Z4.')*p))-(v*((Z4. 1 )*x))); 
phasel_4=exp(-i*k_y*argl_4/c); 

arg2_4=((Y1+Z4).’)*((p-x)+(0.5*(u-v)*(((Yl+Z4).')*(p+x)))) 

phase2_4=exp(i*k_y*arg2_4/c); 

tau_4=l/c*(((Y1+Z4). 1 )*(p-x)); 

neu_4=k_y*(((Y1+Z4). 1 )*(u-v)); 

amb_4=pi*(pulse-abs(tau_4))*exp(- 


i*0.5*neu_4*pulse)*sinc((0.5/pi)*(-neu_4)*(pulse-abs(tau_4))); 
K4(m,n)=phasel_4*phase2_4*amb_4*J; 



% for receiver_5 

argl_5=((Y1+Z5). 1 )*((u*((Z5. 1 )*p))-(v*((Z5.')*x))); 
phasel_5=exp(-i*k_y*argl_5/c) ; 

arg2_5=((Y1+Z5). 1 )*((p-x)+(0.5*(u-v)*(((Y1+Z5). 1 )*(p+x)))); 
phase2_5=exp(i*k_y*arg2_5/c); 
tau_5=l/c*(((Y1+Z5). 1 )*(p-x)); 
neu_5=k_y*(((Y1+Z5). 1 )*(u-v)); 
amb_5=pi*(pulse-abs(tau_5))*exp(- 
i*0.5*neu_5*pulse)*sinc((0.5/pi)*(-neu_5)*(pulse-abs(tau_5))); 

K5(m,n)=phasel_5*phase2_5*amb_5*J; 

% for receiver_6 

argl_6=((Y1+Z6). 1 )*((u*((Z6. 1 )*p))-(v*((Z6. 1 )*x))); 
phasel_6=exp(-i*k_y*argl_6/c); 

arg2_6=((Y1+Z6). 1 )*((p-x)+(0.5*(u-v)*(((Y1+Z6). 1 )*(p+x)))); 
phase2_6=exp(i*k_y*arg2_6/c); 
tau_6=l/c*(((Y1+Z6). 1 )*(p-x)); 
neu_6=k_y*(((Y1+Z6). 1 )*(u-v)); 
amb_6=pi*(pulse-abs(tau_6))*exp(- 
i*0.5*neu_6*pulse)*sinc((0.5/pi)*(-neu_6)*(pulse-abs(tau_6))); 

K6(m,n)=phasel_6*phase2_6*amb_6*J; 

% for receiver_7 

argl_7=((Y1+Z7). 1 )*((u*((Z7. 1 )*p))-(v*((Z7. 1 )*x))); 
phasel_7=exp(-i*k_y*argl_7/c) ; 

arg2_7=((Y1+Z7). 1 )*((p-x)+(0.5*(u-v)*(((Y1+Z7). 1 )*(p+x)))); 
phase2_7=exp(i*k_y*arg2_7/c); 
tau_7=l/c*(((Y1+Z7). 1 )*(p-x)); 
neu_7=k_y*(((Y1+Z7). 1 )*(u-v)); 
amb_7=pi*(pulse-abs(tau_7))*exp(- 
i*0.5*neu_7 *pulse)*sinc((0.5/pi)*(-neu_7)*(pulse-abs(tau_7))); 

K7(m,n)=phasel_7 *phase2_7 *amb_7 *J; 

% for receiver_8 

argl_8=((Y1+Z8). 1 )*((u*((Z8.')*p))-(v*((Z8. 1 )*x))); 
phasel_8=exp(-i*k_y*argl_8/c); 

arg2_8=((Y1+Z8). 1 )*((p-x)+(0.5*(u-v)*(((Y1+Z8).')*(p+x)))); 
phase2_8=exp(i*k_y*arg2_8/c); 
tau_8=l/c*(((Y1+Z8). 1 )*(p-x)); 
neu_8=k_y*(((Y1+Z8). 1 )*(u-v)); 
amb_8=pi*(pulse-abs(tau_8))*exp(- 
i*0.5*neu_8*pulse)*sinc((0.5/pi)*(-neu_8)*(pulse-abs(tau_8))); 

K8(m,n)=phasel_8*phase2_8*amb_8*J; 

K (m, n) = (abs (K1 (m, n) +K2 (m, n) +K3 (m, n) +K4 (m, n) +K5 (m, n) +K6 (m, n) +K7 (m, n) +K8 (m 
,n))).*2; 

end 

end 


Localization in Velocity-Space 

o.g,o,g,o,g,o,g.o.g,o,g,o,g,o,g.o. 

'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O 


nx=length(u_x); 
ny=length(u_y); 
K=zeros(ny,nx); 
Kl=zeros(ny,nx); 
K2=zeros(ny,nx); 
K3=zeros(ny,nx); 
K4=zeros(ny,nx); 
K5=zeros(ny,nx); 


% length of x = columns in K 
% length of y = rows in K 
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K6=zeros(ny,nx); 
K7=zeros(ny,nx); 
K8=zeros(ny,nx); 


for m=l:ny 

for n=l:nx 

u=[u_x_start+((n-l)*res_u);u_y_end-((m-l)*res_u);0]; 

% for receiver_l 

argl_l=((Yl+Zl). 1 )*((u*((Z1. 1 )*p))-(v*((Z1.')*x))); 
phasel_l=exp(-i*k_y*argl_l/c) ; 

arg2_l=((Yl+Zl). 1 )*((p-x)+(0.5*(u-v)*(((Yl+Zl). 1 )*(p+x)))) 
phase2_l=exp(i*k_y*arg2_l/c); 
tau_l=l/c*(((Yl+Zl). 1 )*(p-x)); 
neu_l=k_y*(((Yl+Zl). 1 )*(u-v)); 
amb_l=pi*(pulse-abs(tau_l))*exp(- 
i*0.5*neu_l*pulse)*sinc((0.5/pi)*(-neu_l)*(pulse-abs(tau_l))); 

K1(m,n)=phasel_l*phase2_l*amb_l*J; 

% for receiver_2 

argl_2=((Y1+Z2). 1 )*((u*((Z2. 1 )*p))-(v*((Z2. 1 )*x))); 
phasel_2=exp(-i*k_y*argl_2/c); 

arg2_2=((Yl+Zl). 1 )*((p-x)+(0.5*(u-v)*(((Y1+Z2). 1 )*(p+x)))) 
phase2_2=exp(i*k_y*arg2_2/c); 
tau_2=l/c*(((Y1+Z2). 1 )*(p-x)); 
neu_2=k_y*(((Y1+Z2). 1 )*(u-v)); 
amb_2=pi*(pulse-abs(tau_2))*exp(- 
i*0.5*neu_2 *pulse)*sinc((0.5/pi)*(-neu_2)*(pulse-abs(tau_2))); 

K2(m,n)=phasel_2 *phase2_2 *amb_2 *J; 

% for receiver_3 

argl_3=((Y1+Z3). 1 )*((u*((Z3. 1 )*p))-(v*((Z3.')*x))); 
phasel_3=exp(-i*k_y*argl_3/c) ; 

arg2_3=((Y1+Z3). 1 )*((p-x)+(0.5*(u-v)*(((Y1+Z3). 1 )*(p+x)))) 
phase2_3=exp(i*k_y*arg2_3/c); 
tau_3=l/c*(((Y1+Z3). 1 )*(p-x)); 
neu_3=k_y*(((Y1+Z3). 1 )*(u-v)); 
amb_3=pi*(pulse-abs(tau_3))*exp(- 
i*0.5*neu_3 *pulse)*sinc((0.5/pi)*(-neu_3)*(pulse-abs(tau_3))); 

K3(m,n)=phasel_3 *phase2_3 *amb_3 *J; 

% for receiver_4 

argl_4=((Y1+Z4). 1 )*((u*((Z4.')*p))-(v*((Z4. 1 )*x))); 
phasel_4=exp(-i*k_y*argl_4/c); 

arg2_4=((Y1+Z4). 1 )*((p-x)+(0.5*(u-v)*(((Y1+Z4). 1 )*(p+x)))) 
phase2_4=exp(i*k_y*arg2_4/c); 
tau_4=l/c*(((Y1+Z4). 1 )*(p-x)); 
neu_4=k_y*(((Y1+Z4). 1 )*(u-v)); 
amb_4=pi*(pulse-abs(tau_4))*exp(- 
i*0.5*neu_4*pulse)*sinc((0.5/pi)*(-neu_4)*(pulse-abs(tau_4))); 

K4(m,n)=phasel_4*phase2_4*amb_4*J; 

% for receiver_5 

argl_5=((Y1+Z5).')*((u*((Z5.')*p))-(v*((Z5.')*x))); 
phasel_5=exp(-i*k_y*argl_5/c) ; 

arg2_5=((Y1+Z5). 1 )*((p-x)+(0.5*(u-v)*(((Y1+Z5).')*(p+x)))) 
phase2_5=exp(i*k_y*arg2_5/c); 
tau_5=l/c*(((Y1+Z5).')*(p-x)); 
neu_5=k_y*(((Y1+Z5). 1 )*(u-v)); 
amb_5=pi*(pulse-abs(tau_5))*exp(- 
i*0.5*neu_5*pulse)*sinc((0.5/pi)*(-neu_5)*(pulse-abs(tau_5))); 

K5(m,n)=phasel_5*phase2_5*amb_5*J; 
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% for receiver_6 

argl_6=((Y1+Z6) . ')*((u*((Z6 . ')*p))-(v*((Z6 . 1 )*x))); 
phasel_6=exp(-i*k_y*argl_6/c); 

arg2_6=((Y1+Z6). 1 )*((p-x)+(0.5*(u-v)*(((Y1+Z6). 1 )*(p+x)))); 
phase2_6=exp(i*k_y*arg2_6/c); 
tau_6=l/c*(((Y1+Z6). 1 )*(p-x)); 
neu_6=k_y*(((Y1+Z6). 1 )*(u-v)); 
amb_6=pi*(pulse-abs(tau_6))*exp(- 
i*0.5*neu_6*pulse)*sinc((0.5/pi)*(-neu_6)*(pulse-abs(tau_6))); 

K6(m,n)=phasel_6*phase2_6*amb_6*J; 

% for receiver_7 

argl_7=((Y1+Z7). 1 )*((u*((Z7. 1 )*p))-(v*((Z7. 1 )*x))); 
phasel_7=exp(-i*k_y*argl_7/c) ; 

arg2_7=((Y1+Z7). 1 )*((p-x)+(0.5*(u-v)*(((Y1+Z7). 1 )*(p+x)))); 
phase2_7=exp(i*k_y*arg2_7/c); 
tau_7=l/c*(((Y1+Z7). 1 )*(p-x)); 
neu_7=k_y*(((Y1+Z7). 1 )*(u-v)); 
amb_7=pi*(pulse-abs(tau_7))*exp(- 
i*0.5*neu_7 *pulse)*sinc((0.5/pi)*(-neu_7)*(pulse-abs(tau_7))); 

K7(m,n)=phasel_7 *phase2_7 *amb_7 *J; 

% for receiver_8 

argl_8=((Y1+Z8). 1 )*((u*((Z8. 1 )*p))-(v*((Z8. 1 )*x))); 
phasel_8=exp(-i*k_y*argl_8/c) ; 

arg2_8=((Y1+Z8). 1 )*((p-x)+(0.5*(u-v)*(((Y1+Z8). 1 )*(p+x)))); 
phase2_8=exp(i*k_y*arg2_8/c); 
tau_8=l/c*(((Y1+Z8). 1 )*(p-x)); 
neu_8=k_y*(((Y1+Z8). 1 )*(u-v)); 
amb_8=pi*(pulse-abs(tau_8))*exp(- 
i*0.5*neu_8*pulse)*sinc((0.5/pi)*(-neu_8)*(pulse-abs(tau_8))); 

K8(m,n)=phasel_8*phase2_8*amb_8*J; 

K (m, n) = (abs (K1 (m, n) +K2 (m, n) +K3 (m, n) +K4 (m, n) +K5 (m, n) +K6 (m, n) +K7 (m, n) +K8 (m 

/ n))) . A 2; 

end 

end 


3-D Slice Plot for 

g.g,o,g,o.g.o,g,o,g,o,g,o.g,g.g.o. 

'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o'o 


Position-Space 


p_x_start=-2 0 ; 
p_x_end=2 0 ; 
p_y_start=-2 0 ; 
p_y_end=2 0 ; 
p_z_start=-2 0 ; 
p_z_end=2 0 ; 
res_p=0.5 ; 

p_x=p_x_start:res_p:p_x_end; 
p_y=p_y_start:res_p:p_y_end; 
p_z=p_z_start:res_p:p_z_end; 

nx=length(p_x); 
ny=length(p_y); 
nz=length(p_z); 

K=zeros(ny,nx,nz) ; 

Kl = zeros(ny,nx,nz) ; 

K2=zeros(ny,nx,nz); 


% define vector of x-axis values 
% define vector of y-axis values 


% length of x = columns in K 
% length of y = rows in K 
% length of z = pages in K 
% initialize K matrix 


for a=l:nz 

for m=l:ny 
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for n=l:nx 


p=(10 A 3)*[p_x_start+((n-1)*res_p);p_y_end-((m- 
1)*res_p);p_z_start+((a-1)*res_p)]; 

% for receiver_l 

argl_l=((Yl+Zl). 1 )*((u*((Z1.')*p))-(v*((Z1. 1 )*x))); 
phasel_l=exp(-i*k_y*argl_l/c); 

arg2_l=((Yl+Zl). 1 )*((p-x)+(0.5*(u-v)*(((Yl+Zl).')*(p+x)))); 
phase2_l=exp(i*k_y*arg2_l/c); 
tau_l=l/c*(((Yl+Zl). 1 )*(p-x)); 
neu_l=k_y*(((Yl+Zl). 1 )*(u-v)); 
amb_l=pi*(pulse-abs(tau_l))*exp(- 
i*0.5*neu_l*pulse)*sinc((0.5/pi)*(-neu_l)*(pulse-abs(tau_l))); 

K1(m,n,a)=phasel_l*phase2_l*amb_l*J; 

% for receiver_2 

argl_2=((Y1+Z2).')*((u*((Z2.')*p))-(v*((Z2. 1 )*x))); 
phasel_2=exp(-i*k_y*argl_2/c); 

arg2_2=((Yl+Zl). 1 )*((p-x)+(0.5*(u-v)*(((Y1+Z2).')*(p+x)))); 
phase2_2=exp(i*k_y*arg2_2/c) ; 
tau_2=l/c*(((Y1+Z2). 1 )*(p-x)); 
neu_2=k_y*(((Y1+Z2). 1 )*(u-v)); 
amb_2=pi*(pulse-abs(tau_2))*exp(- 
i*0.5*neu_2 *pulse)*sinc((0.5/pi)*(-neu_2)*(pulse-abs(tau_2))); 

K2(m,n,a)=phasel_2*phase2_2*amb_2*J; 

K (m, n, a) = (abs (K1 (m, n, a) +K2 (m, n, a) ) ) . ^2 ; 


end 

end 

end 


Ambiguity Implementation 

g.g,Q.g,g.g,g.g,o,g,o.g,o,g,o,g,o. 

'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O 


for Chirp Signal 


amb_l=pi*(pulse-abs(tau_l))*exp(i*0.5*(- 
neu_l+(gamma*tau_l))*pulse)*sinc((0.5/pi)*(-neu_l+(gamma*tau_l))*(pulse- 
abs (tau_l))); 


0,Q,O.Q,O.Q,O.Q,O.Q,O.Q,O.Q,O.Q,g. 

'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O'O 

Output Plot 


clear all 

image(p_x,p_y,64*K/max(max(K))) 
hold on 

plot (2,0, 1 +k 1 ,-5,5, 1 A k 1 ,15,15, ,/v k',-2,-15, ,/v k',15,- 
15, 1A k', 1 linewidth 1 , 2 ) 
hold off 

xlabel({ 1 \fontsize{ll) x-component 1 , 1 \fontsize{ll) (position)'}) 
ylabel({ 1 \fontsize{ll) y-component 1 , 1 \fontsize{ll} (position)'}) 
title({ 1 \fontsize{12} Transmitter/Receiver Geometry 1'}) 

% Slice Plot 
norm=max(max(max(K))); 

[UX,UY,UZ]=meshgrid(u_x,u_y,u_z); 
slice(UX,UY,UZ,K/norm, [0], [0], [0]) 
colorbar 
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